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Abstract 

The distance of closest approach of hard particles is a key parameter of their interaction and 
plays an important role in the resulting phase behavior. For non-spherical particles, the distance 
of closest approach depends on orientation, and its calculation is surprisingly difficult. Although 
overlap criteria have been developed for use in computer simulations , no analytic solutions 
have been obtained for the distance of closest approach of ellipsoids in 3-D, or, until now, for 
ellipses in 2-D. We have derived an analytic expression for the distance of closest approach of the 
centers of two arbitrary hard ellipses as function of their orientation relative to the line joining 
their centers. We describe our method for solving this problem, illustrate our result, and discuss 
its usefulness in modeling and simulating systems of anisometric particles such as liquid crystals. 
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I. INTRODUCTION 



Short range repulsive forces between atoms and molecules in soft condensed matter are 
often modeled by an effective hard core, which governs the proximity of neighbors. Since the 
attractive interaction with a few nearest neighbors usually dominates the potential energy, 
the distance of closest approach is a key parameter in statistical descriptions of condensed 
phases. Simple atoms and molecules with spherical symmetry can be viewed as having 
spherical hard cores; the distance of closest approach of the centers of identical hard spheres 
in 3-D or of hard circles in 2-D is the diameter. For non-spherical molecules, such as the 
constituents of liquid crystals, the distance depends on orientation, and its calculation is 
surprisingly difficult [3]. The simplest smooth non-spherical shapes are the ellipse and the 



el 



ipsoid. Although overlap criteria have been developed for use in computer simulations 
, |2j, no analytic solutions for the distance of closest approach have been obtained for 
ellipsoids in 3-D, or, up to now, for ellipses in 2-D. The problem of determining the distance 
of closest approach for two ellipses is particularly intriguing because of its seductive apparent 
simplicity [3j . We have recently succeeded in deriving an analytic expression for the distance 
of closest approach of the centers of two arbitrary hard ellipses as function of their orientation 
relative to the line joining their centers. We describe our method for solving this problem, 
give the solution, illustrate our results, and discuss its usefulness in modeling and simulating 
systems of anisometric particles such as liquid crystals. 



II. STATEMENT OF THE PROBLEM 



We consider two ellipses E\ and E2 in 2D with semi- axes lengths a, and bi where > 
/ b 2 

eccentricity = Wl ^, and major axes oriented along the unit vectors kj (i = 1,2). 

V « 4 

Initially the ellipses are distant so that they have no point in common. One ellipse is then 
translated towards the other along the line joining their centers until they are in point 
contact externally (see Fig. [I]). The problem is to find the distance d between centers when 
the ellipses are so tangent; that is, to find the distance of closest approach. 
The equation of the ellipses are: 

E x : ri A iri = l, Ax = ^1+^-1^!^ =I(I- e 2 kl k 1 ) (1) 
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FIG. 1: Two externally tangent ellipses E\ and E%. The directions of the major axes are given by 
ki and k 2 , the centers are joined by the vector d. 

and 

E 2 : r 2 A 2 r 2 = 1, A 2 = 1 ^1 + ^ - 1^ k 2 k 2 ^ = i(I - e 2 2 k 2 k 2 ), (2) 

where I is the identity matrix and kjkj is the dyad product. The vector joining the centers 
is given by d = dd; d is a given unit vector. Our goal is to find the distance d as function 
of ellipse parameters ai, fei, a 2 , b 2 and orientations ki • d, k 2 • d and ki • k 2 . 
It is tempting to seek a solution by solving the quadratic equations 

riA iri = 1 (3) 

and 

(n - d)A 2 (n - d) = l (4) 

simultaneously for the points of intersection, and then requiring that the distance d between 
centers be such that there is intersection exactly at one point. This approach fails for the 
following reason: although the components of ri at the points of intersection can be obtained 
by solving a quartic equation (say for the x-component of ri), the condition requiring that 
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the quartic have exactly one double real root is not straightforward to implement (there 
are four roots, and it is not clear which two roots need to coalesce to yield the required 
tangency condition) and it further gives an equation in d whose order is higher than quartic, 
and which cannot therefore be solved analytically. 



III. THE SOLUTION 



Our approach proceeds via three steps: 

(1) Transformation of the two tangent ellipses E\ and E 2 , whose centers are joined by 
the vector d, into a circle C[ and an ellipse E 2 , whose centers are joined by the vector d'. 
The circle C[ and the ellipse E' 2 remain tangent after the transformation. 

(2) Determination of the distance d! of closest approach of C[ and E' 2 analytically. 

(3) Determination of the distance d of closest approach of E\ and E 2 by inverse transfor- 
mation of the vector d'. 



A. Transformations 




FIG. 2: The transformed circle C[ and ellipse E 2 after applying transformation T to the ellipses 
in Fig. [TJ The major axis of transformed ellipse E 2 is along k , and the centers are joined by the 
vector d'. 
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An ellipse can be transformed into a unit circle by anisotropic scaling. We introduce for 
this purpose the matrix T, which transforms the ellipse E\ into a unit circle C[ and the 
ellipse E 2 into another ellipse E 2 . The transformation is a scaling by the factor l/a\ along 
the kx direction and by the factor 1/bi in the direction perpendicular to ki. The transfor- 
mation matrix T, which transforms position r to a position r' in a space with dimensionless 
coordinates, is 

and the inverse, T _1 , is 

T-^&i^ + T&kx), (6) 

where 



i=lr-i)- P) 



One can easily verify that T _1 AiT _1 = I. 

If = Tr-j, or, equivalently, = T _1 r^ (i = 1,2), substitution into Eq. (JTJ) gives a unit 
circle C( and a new ellipse I<7 2 (see Fig. [2]). That is, 

: riAiri = riT^Ail- 1 ^ = = 1 (8) 



and 



A' can be written as 



E' 2 : r 2 A 2 r 2 = ^T~ 1 A 2 T~ 1 i^ = r' 2 AY 2 = 1. (9) 



A' = |(I + r/kxki)^ - e^k 2 k 2 )(I + r/kxki). (10) 

The eigenvectors of A' provide information about the directions of the principal axes and 
the eigenvalues about the lengths of the semi-axes of the transformed ellipse E' 2 . Since A' 
is real symmetric, its eigenvalues A + and A_ are real, and the corresponding eigenvectors 
k^_ and k_ are orthogonal. Explicit expressions for these are given in the Appendix. The 
lengths of the semi-axes, if A + > A_, are given by 

a' 2 = 

and we note that a! 2 > b' 2 . 



A_i_ 

(11) 
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Under the transformation T, the vector d is transformed to 



d' = Td = dTd = d'd r , (12) 



where d = — — is a unit vector. Explicitly, 
lTdl 



and 



Td^Ifd+f^-ll^-d)^) (13) 



|Td| = l^/l-efCki-d) 2 . (14) 
B. Distance d' of closest approach of a circle and an ellipse 

We next derive the a useful relation between the position vector r of a point on the ellipse 
and the unit outward normal n at that point. For an ellipse, given by rAr = 1, the unit 
normal n is 

„ _ V(rAr) Ar 

n ~ |V(rAr)| ~ v^A^' (15) 
Multiplying Eq. (03J) by 1 = A" 1 gives 

Bn = J- : (16) 

VrA 2 r 

and multiplying Bn by n gives 

nln = — — . (17) 
rA 2 r v ; 

Substituting into (flBT) . we obtain r in terms of the unit normal n 



(18) 



VnBn 

If a unit circle and an ellipse are externally tangent, then the directions of their normals 
at the point of contact must be opposite. If the unit outward normal of the unit circle C[ 
at the point of contact is n', then 



r\ = n', rf, = (19) 



and we have, for the vector d' joining the centers, 



where 



/2 



6 2 ( I + 5k'_k'_ ) (21) 

and 

5 = ^ - 1 > 0. (22) 

2 

Eq. f[2"Uj) is a key result. It is a vector equation with only two unknowns: the magnitude of 
d' and the direction of n'. It can be solved for d! as follows. We multiply both sides of Eq. 
( 120]) by k'_ and by k' + , and letting 

k_ ■ d = sin0, k' ■ n' = sin?/), k' + ■ d = cos</>, k' + ■ n' = cos?/', (23) 

we get, from Eq. (1201) . 

d'sin0 = sinV | 1 + —, b ML=L= ) (24) 
\ Vl + 5 sin 2 ^ J 

and 

d'cosift = cos?/; [ 1 + — = &2 1 . (25) 

Here the unknowns are ip and d! . In the special case of 5 = 0, d! = 1 + 6 2 = 1 + a' 2 , and in 
the case of = — , d! = 1 + b 2 \fL + 5 = 1 + a' 2 . In general, ^ — , and the solution for <f is 
more challenging. 

We let q = \/T+~6 sin 2 ip, then 

sin 2 jj, = <tzl (26) 
o 

and 

a 2 - I 

cos 2 ?/' = 1 - . (27) 

o 

Substitution into Eq. (1241) (a-b). squaring both sides and dividing these two equations to 
eliminate d! gives a quartic equation for q, 

tan 2 <f>(5 + 1 - q 2 )(ji + l) 2 = (q 2 - 1)(£ + 1 + 5) 2 . (28) 

°2 °2 

This can be written in the standard form Aq 4 +Bq 3 +Cq 2 +Dq+E = 0, where the coefficients 
are 

A = — ^(1 + tan 2 0), (29a) 

^2 

B = -^-(1 + tan 2 + 5), (29b) 



C = - tan 2 <f> - (1 + S) 2 + \{l + (1 + S) tan 2 0), (29c) 
D = |-(l+tan 2 0)(l + 5), (29d) 

"2 

E = (l + tan 2 + 5)(l + 5) (29e) 

and 

(k' -d') 2 , , 

tan 2 = — ^ 4—. 30 

l-(k'_-d') 2 

The roots of Eq. (125)1 can be obtained explicitly as follows. 

To make contact with the standard solution of the quartic equation, using Ferrari's 
method |4|, we define 

3B 2 C , . 

a = -8Ai + A> (31a) 

-35 4 CB 2 BD E 



and 



and 



^ = -^"7, (31d) 

Q = -^ + ^-^ (31e) 
^ 108 3 8 v ; 

where we take the principal values of the roots. If U — 0, then 

y = ~a-0 1/S , (32) 

otherwise 

In terms of these, the one real positive root q is 



q = -^ + U^T2-y + J-(3a + 2y +7 ^ Tv )\. (34) 
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In the special case when a + 2y = 0, then (3 = (which we have not observed in this 
problem, but include here for completeness) and the positive real root is given by 



q = -AA + i 2 • (35) 

Knowing q, d! can be found by squaring both sides of Eqs. f l24l) (a-b) and adding; this 
gives 



d , ^ l+ m±m\ u_t=i (36) 



2 

8 V Q J V 5 J V ' Q 

The vector joining the centers of the circle and the ellipse is given by 



d' = d'd'. (37) 



C. Distance d of closest approach 

The distance of closest approach of the two ellipses is obtained via the transformation 
from d' to d, 

d = T-M' = d'T-'d' = d'T- 1 ^- = -id = - — d ' H (38) 

ITdl lTdl 1 . A 7J7 ^2 



-Vl-e?(k!-d) 
and finally we have 

d ' h. (39) 

1-eKkx-d) 2 

This is the solution for the distance of closest approach, which is our main result. 



D. Contact point 

In addition to the distance of closest approach, it is interesting and useful to locate the 
point contact. We denote the vector from the center of Ellipse 1 to the point of contact as 
r c . In the transformed coordinate system, where Ellipse 1 has become a unit circle after the 
affine transformation, the vector from the center of circle to the point of contact is r' c = n' 
where n' is the unit normal at the point of contact, whose components along the orthogonal 
unit vectors and k'_ are known Thus, to obtain r c , it is only necessary to perform the 
transformation of n', that is, 

r c = J- l n (40) 
9 



The components of n' are given by 



k'_ • n' = sin-0, k' + • fi' = cos^ 

where ip is given by 



sinip = sgn(sin(f))\/ — - — , 



q 2 - 1 

cos^ = sgn(cos(f)) d 1 — . 

where sgn(x) gives the sign of x and the angle is known. 
It follows that 

n' = cos + sin ipk.'_ 
Writing k' + and k'_ in terms of k\ and k 2 

f/ (ki + k 2 ) , . (ki-k 2 
k, = cos 7 = + sin 7 



where 



and 



2 v / l + k 1 -k 2 y/2y/l - ki • k 2 

ki + k 2 ) 



cos 7 = k'. 



sin 7 = k. 



V2Vl + ki • k 2 
(ki - k 2 ) 



2Vl-ki-k 2 



f , . (ki + k 2 ) (ki-k 2 ) 
k = — sin 7 = + cos 7- 



V2Vl + ki • k 2 V2Vl - ki ■ k 2 

Substitution gives n' in terms of k x and k 2 gives 



-/ / / , \ (ki + k 2 ) . (ki - k 2 ) 

n = cos(^ + 7) == =f= + sm(i/) + 7)- 



y/2 y/l + kx • k 2 v^V 7 ! - ki • k 2 

Now 

r c = T- 1 n' = 6i (l + 77kiki) n' 

and so 

7 / / , v (ki + k 2 ) . (ki - k 2 ) 

r c = & lC os(^ + 7) ~ j 5 == + &i sm(^ + 7) ~ , 5 5 + 
v^V 1 + ki • k 2 V2V 1 - ki • k 2 

(l + ki-k 2 )ki . (1 — k 

(ai - 61) cos(^ + 7) ~ ; 5 ; + (ai - &i) sm(^ + 7) - - 

V2V1 + k x • k 2 v^V 1 
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and finally 

u , , s," <■ \ cos(-?/> + 7) , , , nT r n sin(^ + 7) f 
r c = [( ai + (ax - 6i)k a ■ k 2 ) - AI_^= 4= + (ai - K - 61)^ • k 2 ) - ^= ^4=]^ + 

v^V 1 + ki • k 2 V 1 - ki • k 2 

1 1 

[6icos(^ + 7) . . -6isin(^ + 7) 5 . ]k 2 (52) 

V^Vl + kx • k 2 V2V1 - ki • k 2 

IV. DISCUSSION 

The above methods give a closed form expression for the distance of closest approach 
and the position of the point of contact for two ellipses of arbitrary size, eccentricity and 
orientation. Detailed steps of the calculation are given in the Appendix. To demonstrate 
the applicability of the method, we give two examples: calculation of the excluded area and 
the locus of the point of contact while one ellipse is fixed and the other is rotated. 

A. Excluded area A ex 




FIG. 3: Excluded area for two ellipses. Ellipse E\ is fixed at origin, and ellipse Ei rotates around 
it, keeping its orientation fixed and remaining tangent to E\. The center of E% traces out the 
dashed curve. The area bounded by the dashed curve is the excluded area A ex . 

From the analytical solution provided in Section 1111} one can easily compute, numeri- 
cally, the excluded area for two identical ellipses whose orientation is fixed by integrating 
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d 2 (d, ai, 61, a 2 , 62, ki, k 2 ) over d; 

4x = ^ y £Z 2 (d,ai,fe 1 ,a2,fe2,ki,k 2 )|dd|. (53) 

Fig. [3] shows the locus of the center of ellipse E 2 rotating around E\ while keeping the 
orientation of both ellipses fixed. Here a\ — a% = 2, b\ = 62 = 1- When the angle between 
the major axes is 30°, the excluded area is 26.4 (Fig. Oa). If the angle is increased to 45°, 
then the excluded area is 27.6. If the angle is 90°, then the excluded area is 29.7 (Fig. 
Ob). The excluded area increases monotonically with the angle between major axes of two 
ellipses; it is the smallest when the major axes are parallel, and the largest when the major 
axes are normal to each other. 



B. Locus of the point of contact 




FIG. 4: Locus of the point of contact. Ellipse E\ is rotated about its center, while ellipse E2 keeps 
its orientation fixed. The center of E2 moves so that Ei remains tangent to E\. The point of 
contact traces out the dashed curve. 

Fig. H] shows that locus of the point of contact when ellipse E\ is rotating about its center 
while ellipse E 2 keeps its orientation. It is interesting and unexpected that the locus has 
dipolar rather than quadrupolar symmetry. 
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C. Potential Applications 



Our result, the analytical expression for the distance of closest approach of two hard 
ellipses, has a number of potential applications. It may be useful in modeling 2-D liquid 
crystals, both analytically and numerically The excluded area, discussed above, is a key 
parameter in statistical models 5j which can be calculated for ellipses from our result. 
Another potential application is in the theory of nematic liquid crystals. One important 
contribution to the elastic constants of nematics is due to anisotropic dispersion forces. The 
average Van der Waals interaction energy of a molecule with its neighbors is an algebraic 
function of the distance of closest approach. The origin of three distinct elastic constants in 
nematics is still unresolved. Our result may be useful in modeling elastic constants in 2-D 
nematics, and possibly giving insights towards understanding their origins in general. 

Monte Carlo calculations have played an important role in modeling the phase behavior 
of isotropic fluids and liquid crystals [6( . Vieillard-Baron developed the first overlap criterion 
for identical hard ellipses He derived a contact function \l/(a, b, k\, k2, d) such that ^ = 
when the ellipses are tangent (either exteriorly or interiorly), and this function is positive 
and at least one of two auxiliary functions are negative if the ellipses have no real point in 



common. 



systems 



EL 



Tris overlap criterion has been used in Monte-Carlo simulations of hard ellipse 



, It may be possible to solve ^(a, b, k 1; k 2 , d) = for d (this involves solving a 
quartic equation), and thus obtain a result similar to ours; to our knowledge this has not yet 
been done. However, Vieillard-Baron's contact function \1/ is only valid for identical ellipses, 
and so this result would not be as general as ours, presented here. 

According to the Hohenberg-Mermin- Wagner theorem, long range order corresponding 
to broken continuous symmetry is not allowed in 2D systems with short-range interactions 
js]. The possibility of long-range order in 2D nematics has been discussed theoretically [s| 
and examined using Monte- Carlo simulations with Lennard- Jones like potentials [lOj] as well 
as with hard rods Uj . Although it has been shown that true long range order cannot exist 
if the interparticle potential is separable into a positional and an orientational part 10], it 
is not clear what the implications are for systems of hard ellipses. Frenkel has shown that 
only quasi-long range order exists for hard spherocylinders ll|, that is, the correlations 
in orientational order decay algebraically. Hard ellipsoids, however, can show dramatically 
different behavior from hard spherocylinders 12| (hard ellipses do not form smectic phases, 
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whereas spherocylinders do [13{]), and for this reason Monte Carlo simulations of hard ellipses, 
on systems larger than studied by Vieillard-Baron, would be of considerable interest. Our 



result for the distance provides an overlap criterion which could be useful 



y applied here. 



Another area of interest is phase separation in hard particle systems |14l . Ilq . For example, 
simulations of hard disks and hard parallel squares have been studied, and phase separation 
has been observed. Theoretical studies, on the other hand predict no phase separation in 



2D 



161 ] . Our results could provide the criterion for the overlap of ellipses of different sizes, 



and thus enable Monte-Carlo simulations of binary mixtures of hard ellipses. 

Vieillard-Baron also provides an overlap criterion for two identical ellipsoids of revolution 
in 3D [l|. This involves the evaluation of a contact function \l/ and five auxiliary functions, 
three of which must be non-negative and at least one among the remaining three must be 
negative to avoid overlap. Perram and Wertheim provided a more general overlap criterion 
for hard ellipsoids [2]. Their scheme for evaluating the criterion involves an iterative nu- 
merical technique to find the maximum of a scalar function. Our results can provide the 
basis of a simple algorithm to determine the distance of closest approach of two ellipsoids in 
3D. This involves passing a plane through the line joining the centers of the two ellipsoids, 
determining the distance of closest approach of the ellipses in the plane, then rotating the 
plane and finding the largest such distance. The details of this algorithm will be published 
elsewhere. 



V. CONCLUSION 



We have derived an analytic expression for the distance of closest approach of two hard 
ellipses with arbitrary orientation in 2D. The strategy is to transform the ellipses into a 
circle and a new ellipse by a scaling transformation. The relation between the position of 
a point on the ellipse and the normal at that point allows the tangency condition between 
the circle and ellipse to be written as a simple vector equation with two unknowns, which 
may be solved analytically for the distance between the centers. The solution requires the 
solution of a quartic equation, whose single positive real root can be uniquely determined. 
The final result for the distance is obtained by the inverse scaling transformation. Explicit 
instructions for calculating the distance are given in the Appendix. Our result may be useful 
in analytic and numerical models of orientationally ordered systems. 
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VI. APPENDIX 



A. Detailed calculation of the distance of closest approach 

We start by specifying all the quantities required for the calculation of d. 
ai and bi are the lengths (a\ > b\) of the major and minor axes of ellipse E\. 
a 2 and b 2 are the lengths (a 2 > b 2 ) of the major and minor axes of ellipse E 2 . 

ei — \ jl 1 is the eccentricity of ellipse E 1 . 



e 2 = \ / 1 1 is the eccentricity of ellipse E 2 . 

V °2 

k x • d is the cosine of the angle between the major axis k x of ellipse E x and the direction 
d of the line joining the centers. 

k 2 ■ d is the cosine of the angle between the major axis k 2 of ellipse E 2 and the direction 
d of the line joining the centers. 

k x ■ k 2 is the cosine of the angle between the major axis k\ of ellipse E 1 and the major 
axis k 2 of ellipse E 2 . 

The above quantities are specified in the statement of the problem. 

The following quantities are derived from these. 

V = y i -1- (54) 

In the coordinate system with the basis (ki +k 2 ) /y/2~+ 2k\ • k 2 and (ki — k 2 )/y/2 — 2ki ■ k 2 , 
the components of A' are 

A' n = |(i + \{i + ^ ■ k 2 )(^( 2 + V )- 4(1 + • k 2 )) 2 )), (55) 



A ' 22 = l (1 + i (1 " fcl ' ^ 2){r]{2 + r]) ~ 62(1 " ^ ' fc2))2)) (56) 



1 



and 



bjl 



A' 12 = A 21 = ^- yj 1 - (kx • k 2 )*fa(2 + rj) + e 2 (l - rffa • k 2 ) 2 )). (57) 
The eigenvalues of A', in terms of these, are 

A+ = \(Ki + A' 22 ) + ^\(A' n - A' 22 y + A 



,/2 
H2 



(58) 



15 



and 



It follows that 



A- = \(Ki + A' 22 ) - J\(A' U - A' 22 f + A; 2 2 . (59) 



b' 2 = ^= (60) 

A4 



a' 2 = * (61) 

A_ 



and 



The eigenvectors are given by 

^ = 1 ( A; 2 (k 1 + k 2 ) + (A + -A; 1 )(k 1 -k 2 ) (g2) 



and 



Then 



V2\Ja[ 2 2 + (A+ - A' n ) 2 V V 7 ! + ki • k 2 V 7 ! - ki • k 2 

(A+-A' 11 )(k 1 + k 2 ) , A' 12 (k!-k 2 ) 



V2\/A' 12 2 + (A+ - A' n ) 2 V V 7 ! + ki • k 2 v 7 ! - k x • k 2 



k' • d' = cos = x 



v^A;, 2 + (A + - A' n y yj 1 - e?(ki • d7 



kx • d) + (k 2 • d)+ - - 1 (kx • d)(kx • k 2 ))+ 



v'l • k, • k, "■■ ' ' \ a i 

(A+-A' n ) f b lf f 3^ n 3 \ f&i 



(63) 



(^(kx • d) - (k 2 • d)- (*l - 1) (kx • d)(kx • k 2 )) ] . (64) 

ax \ a i / / 



V 7 ! - ki • k 2 ai 

If kx = — k 2 , then — k 2 may be replaced by +k 2 without the loss of generality. If kx = k 2 , 
care must be taken evaluating the above expression. Letting kx • k 2 = cos 9, in the limit as 
9 -> we find that if A' n > A 22 , then (A+ - A' n ) ~ Af 2 , and 

k' + -d' = cos0 = 1 — (ki-d), (65) 

Vl-e^kx-d) 201 

while if A' u < A' 22 , then 



'i-(kx-d) 2 

k + • d — cos <p — — ■■ . (66) 



l-ef(kx-d) 
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Next, 

(67) 

If = it/ 2 or if 5 = 0, then cf = 1 + a' 2 , and Eq. (!73l) can be evaluated directly. Otherwise, 

tan 2 = -s— — s 1, (68) 

(k' + -d') 2 

A = — ^(l + tan 2 0), (69a) 

B = -^-(1 + tan 2 + 5), (69b) 

C = _ tan 2 - (1 + 5f + (1 + (1 + 5) tan 2 0), (69c) 

D= ^-(l+tan 2 0)(l + <5), (69d) 
£ = (1 + tan 2 + 5)(l + 5), (69e) 



and 



If ^ 0, then 



a = "8^ + A (7 ° a) 



„ fi 3 BC D , . 



_ -3ff 4 Cff 2 _BD_ E_ 
7 ~ 256I 1 + 16l3 ~ 4^2 + A' { ] 

^ = -^-7, (70d) 

Q = -^ + ^~^ (70e) 

1/3 



a 3 


«7 


/? 2 


108 


8 ' 


i+ 




27 J 



and the principal values of the roots are taken throughout; 

if U = 0, 

6 
and 




(70s 



^ + i ( \/a + 2j + J- (3a + 2y+ -y=5j ) I- (70h) 
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If (3 = 0, then 

" = " 1.1 



q= b + .=a±^zi2, (71 , 



and finally 



61. (73) 



^l-e?(ki-d) 

B. Detailed calculation of the position of the point of contact 

To obtain the vector r c from the center of Ellipse 1 to the point of contact, we first need 
to compute n' and then perform the inverse afline transformation. Explicitly, r c = T _1 n'. 
The components of n' can be defined though the inner product with k' and k'_:: 



q 2 — 1 

k' + ■ n' = cos^ = sgn(cos(f))d 1 — , (74) 

- 1 

k_ ■ n' = sin^ = sgn( sin 0)y — - — . (75) 

where sgn(x) gives the sign of x. The expression for cos0 is given in Eq. ??. sin0 can be 
calculated similarly, and 



k' • d' = sin = x 

v^7 A i2 2 + (A + -A' n )V 1 - eKk ■ df 

(A+ " Ail) ' W % ■ d)+(k 2 • d)+ - l) ( kl ■ d)( kl ■ k 2 ))+ 



\fl + ki • k 2 °i " V«i 

A' 12 /^i/r i\ /r In /^i 



fix • d) - (k 2 • d)- (h. - 1^ ( kl ■ d)( kl ■ k 2 )) j . (76) 



\A - ki • k 2 °i 
It follows that 

n' = cos ^k' + + sin ^k'_ . (77) 
Writing k' + and k'_ in terms of ki and k 2 

f/ (ki + k 2 ) . (ki - k 2 ) 

k + = cos 7 = + sin 7 =, (78) 

v 7 ^ V 1 + ki • k 2 V2 v 1 - ki • k 2 



18 



t , . (kl + k 2 ) (kl - kg) 

k_ = -sm7 = — + cos7 = =, (79) 

>/2 V 1 + k x • k 2 V2 V 1 - ki • k 2 



where 



and 



oos7 = kV- / £l+fe2) . , (80) 
v^V 1 + ki • k 2 

sin 7 = k- + . . (81) 

V2 VI - ki • k 2 

The angle 7 can be calculated explicitly, since k' + and k'_ are known. Substitution gives 
n' in terms of ki and k 2 and gives 

n' = cos(^ + 7) (kl + i + sin(^ + 7)—^^==- (82) 
V2V 1 + ki • k 2 V2V 1 - ki • k 2 

Finally, after the transformaton T _1 , the contact point r c is given by a linear combination 
of ki and k 2 : 



r c = T _1 n' = 61 (i + r^kiki j n' 



u 1 r \,~ r \ cos(-^ + 7) , , . x sin(-^ + 7) ,, 

= [(ai + (ai - 6i)ki • k 2 ) - ^= ^4= + (ai - (ai - 6i)ki • k 2 ) - A_^= 4=]ki + 



[61 cos(-^ + 7)- 



2Vl + ki -k 2 

== - bi sin('0 + 7)- 



2V1 - ki k 2 
k 2 . (83) 



V2 v 7 ! + ki • k 2 ' V2 V 7 ! - ki • k 2 

There are some special cases, when either q is not given or <p is not well defined We treat 



these separately below. 



I. If S = or 



7r/2, then cos^ and sin^ are not required, since that n' = d'and 
T _1 Td d 



T -i d / 



Td 



^l-e^ki-d) 2 



(84) 



2. If ki = k 2 , care must to be taken in evaluating the angle 0. Since ki = k 2 , these 
vectors do not span the space, and a new vector kf, perpendicular to ki, needs to be 
introduced. . 



19 



(a) If An > A 2 2, then k' + = k l5 and k' = k^, and cos0 and sin0 are given by 

k' + • d' = cos0 = 1 = — (ki • d), (85) 

Vl-eK^-d) 201 

k'_ d' = sin0 = 1 = =(kj- • d), (86) 

Vl-ei(ki-d) 2 

Here n' = cos^kx + sin^kj 1 , and the r c is given by 

r c = T _1 d' = &i (i + r/kiki ) (cos^ki + sm^kf) 

= ai cos^ki + b\ sin^kj 1 . (87) 

(b) If An < A 22 , k' + = kj;, and k'_ = k\, then 

k' + • d' = cos0 = 1 = =(kj L • d), (88) 

Vi- e ?(£i-a) 2 

k'_ • d' = sin0 = 1 = — (ki • d), (89) 

Vl-e^^-d) 201 

and n' = cos ^Iq- + sin ^ki We then have 

r c = T _1 d' = 61 (l + 77k! k^ (cos^kf + sin^kx) 

= bi cos-^kj 1 + ai sin-^kx. (90) 
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